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SUMMARY 

The  Report  describes  a  way  of  designing  a  two  dimensional  array  to 
minimise  the  side-lobes  while  maintaining  any  givnn  beamwidth  of  the  radiation 
pattern  in  the  array  plane.  The  amplitudes  and  phases  associated  with,  the 
elements  are  adjusted  by  a  modified  gradient  method  which  uses  a  linear  pro¬ 
gramming  procedure.  An  example  is  given  in  which  the  side-lobe  level  for  a 
six-element  array  is  lowered  by  21  dB. 
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1  INTRODUCTION 

.  1.1  We  are  concerned  with  arrays  consisting  of  a  fairly  small  number  (K)  of 

elements  placed  in  a  horizontal  plane.  Each  element  is  omnidirectional  and 
4  unaffected  by  its  neighbours.  The  array  polar  diagram  will  be  considered  in 

the  horizontal  plane  only,  and  at  a  single  frequency. 

Having  fixed,  arbitrarily,  the  positions  of  the  elements,  then  for  i 
linear  method  of  beam-forming,  we  have  to  choose  2N  quantities,  namely,  ••• 
amplitude  weighting  and  phase  shift  to  be  applied  to  each  element.  We  noiraPv 
desire  that  the  final  polar  diagram  should  have  a  narrow  main  beam  and  low 
side-lobes.  The  present  work  deals  with  the  problem  which  may  be  stated 
formally  as  follows: 

O' 

For  a  given  beamwidth,  choose  the  weights  and  phases  to  ni i r.  fi  .i  .*;• ■ 
largest  side-lobe.  \ 

1.2  The  Dolph-Tschebyschef f  theory  solves  this  problem  for  linear  and 
equally  spaced  arrays.  However,  as  far  as  I  know,  there  is  no  theory  available 
for  the  present  case.  The  brute  force  method  of  trying  all  combinations  of 
weights  and  phases  is  quite  impossible,  unless  N  is  trivially  small,  since*  it 
oquires  far  too  much  computer  time.  The  only  practicable  method  seems  to  bo 
to  select  an  arbitrary  initial  situation  and  then  try  to  reduce  the  sio*  lol-.* 
love!  by  some  kind  of  gradient  method. 

In  section  2  the  method  used  in  this  Report  is  explained.  It  will  be 
seen  to  be  a  gradient  method,  but  with  special  treatment  to  deal  with  discon¬ 
tinuities  in  the  gradient  of  the  object  function.  The  special  treatment  a 

linear  programming  procedure.  In  section  3  the  method  is  applied  to  -u.-  -<y 

under  consideration.  In  section  4  numerical  results  for  a  particular  array  are 
presented . 

2  GENERAL  METHOD 

2.1  In  this  section  we  present  in  general  terms  the  method  of  solution. 

Any  set  of  weights  and  phases  may  be  represented  as  a  vector  in  a  space 
of  2N  dimensions.  Since  the  array  pattern  should  have  a  preset  beamwidth, 

•  the  vectors  are  restricted  to  a  space  of  dimension  n,  say,  which  is  loss  than 

2N.  It  might  be  expected  that  we  could  choose  n  coordinates  arbitrarily,  and 
,  solve  for  the  remaining  2N  -  n.  This  tur.-.s  out  to  be  the  case,  ;>o  we  can 

consider  the  problem  to  be  one  of  minimising  the  side-lobe  level  when  the  point 
of  interest  is  allowed  to  range  freely  over  a  space  of  n  dimensions. 
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2.2  We  shall  denote  the  point  of  interest  by  z  ■  (z^,  ...  ,  z^) ,  the  number 
of  side-lobes  by  tn  and  the  values  of  the  power  at  these  side-lobes  by 

^l’  *  ^m*  these  are  functions  of  £.  It  may  be  that  for  some  z, 

m  ■  0.  There  are  then  no  side-lobes  and  the  problem  is  solved.  In  general 
however  m  >  0  and  so  there  exsts  a  largest  y^  with  value  y: 

y  -  max  (y^)  .  (1) 

i  * 

In  order  to  reduce  y  it  is  natural  to  make  a  step  in  the  direction  oi 
-grad  y 

Az.  ■  -  c  grad  y  (2) 

where  c  >  0.  The  simplest  rule  is  to  set  e  equal  to  a  constant,  and  in  fact 
the  program  to  be  described  uses  just  this  method  in  its  initial  stages.  The 
flow  chart  is  shown  in  Fig.i. 

2.3  After  a  number  of  cycles  of  this  simple  gradient  method,  the  procedure 
gets  into  difficulties.  This  happens  when  two  or  more  side-lobes  are  nearly 
equal,  and  the  identity  of  the  largest  changes  from  cycle  to  cycle.  The 
gradient  vector  of  y  is  discontinuous.  It  can  happen  that  y  is  actually 
increased  by  a  gradient  step,  due  to  the  role  of  maximum  side-lobe  passing  from 
one  lobe  to  another.  In  the  program  it  was  arranged  that  if  y  failed  to 
decrease  over  three  successive  iterations,  the  simple  gradient  method  would  be 
abandoned  in  favour  of  the  linear  programming  method  described  below. 

2.4  In  the  neighbourhood  of  the  point  of  interest  ,  we  may  suppose  that 
the  following  linear  approximations  hold 

y.(z)  -  y.(z  )  ♦  (z  -  z  )  .  (grad  y.)  .  (3) 

1  —  1  “O  —  — "O  1  L  r\ 

Thus  the  change  in  y^  is  proportional  to  the  projection  of  the  displacement 

(z  -  z^)  on  the  gradient  grad  y^.  In  order  to  change  all  the  y^s,  it  is 

most  efficient  to  use  a  displacement  which  belongs  to  the  m-dimensional  linear 

manifold  spanned  by  the  gradients,  that  is,  fo.  some  scalars  c. ,  ...  c  , 

'  m 
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or 


(*  -  *  ). 


m 

y 

u 

j-1 


c. (grad  y.)j 


(A) 


z  -  z 
—  — o 


m 

y 

j-i 


c.  grad  y. 
i  l 


(5) 


For,  given  any  displacement,  wo  could  obtain  the  same  changes  in  the  y^s  hy 
using  a  smaller  displacement,  namely,  that  obtained  by  projecting  the  given 
displacement  on  to  the  manifold  spanned  by  the  gradients. 


The  point  is  that  instead  of  dealing  with  an  n-dimensiona'  'isplacement , 
we  need  only  consider  the  m-dimensional  vector  £.  Normally  m  is  less  than 
n,  and  so  the  amount  of  computation  will  be  reduced.  In  matrix  notation  we 
put 


z  -  z  »  £  (6) 

where  £  is  an  m  *  n  matrix  whose  ith  row  is  grad  y^  written  as  a  row 
vector,  and  c  is  a  column  vector  of  order  m.  The  original  linearised 
equation  may  be  written 


1  -  k  ♦  &  <*  -  (7) 

where  ^  and  b  are  column  n-vectors  consisting  of  y^(£)  and  y^r^). 
Substituting,  we  obtain 


2  -  +  h  £  (8) 

where  h  is  the  m  x  m  matrix  £  £' . 

In  order  to  avoid  trouble  with  non-linearity,  let  us  restrain  c  so  that 
none  of  its  components  can  numerically  exceed  a  given  quantity  e;  that  is, 

|ci|  <  c  ,  i  «  1,2,. . . ,m  .  (9) 


These  can  be  converted  into  one-sided  constraints  by  setting 
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ci  “  pi  "  qi  1  i  -  1 . m  d°) 

where  the  p. ,  q^  satisfy 

0  *=■  p .  ^  >  ,  0  q  ^  t  ;  i  =  l , .  , .  ,  n>  .  (11) 

Further,  we  introduce  r.,  s.  sueli  that 
’  t  i 

p.  +  r^  ■  e  ,  q^  +  =  e  ;  i  «  1 , . . ,  ,m  (12) 

where  r.^O.s.^Oj  i  =  1 , , , .  ,m. 

i  i  1  1 

We  now  introduce  yet  more  variables  x..  ...  ,  x  such  that 

l  m 

y  “  yi  "*■  xi  *  i  -  1 , .  .  .  ,m  (13) 

Since  y  »  max  (v.),  then  x.  0  for  each  i.  It  is  not  convenient  to  have 

-  i  1  1 

y  appearing  in  more  than  one  equation,  so  we  shall  retain  the  equation 
obtained  from  i  ■  1,  namely 

y  ■  y,  ♦  Xj  (H) 

and  eliminate  y  from  the  rest 


0  = 


>’i  •  >1 


+  x .  -  x. 


(15) 


Substituting  for  y^  the  value 


b .  ♦  y  h. .  c .  (16) 

1  j-1  1J  J 


as  given  by  the  matrix  equation  (8),  we  obtain 
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ill 

y  ■  b,  +  l  h, .  c.  +  x. 

1  lj  J  1 


O  -  b.  -  b, 

i  1 


in 

*  *i  ■  xi  *  £  <hij  -  v  cj  • 


i  =  2  ,  . .  .  ,m  (1  8) 


y  ■  bi  *  I  hLi  ^  +  xi 

j»i  j 


0  -  b.  -  b,  +  l  (hj,  -  h..)  (p.  -  q.)  +  x.  -  x  ;  i  -  2,...,in(20) 

j«l  J  J  “  * 


2.5  To  summarise,  we  now  have  a  standard  problem  in  linear  programming, 
involving  the  (2m  +  1)  basic  variables 


X1  1  P1 .  Hm  ’  Hl*  ••• 


Pm  »  'll  * 


and  the  (3m  -  1)  alack  variables 


2’  *  *  *  ’  m  ’  ‘1’  ’  m  ’  °1 


r  )  s.  i  « •  •  )  s 


These  are  linked  by  the  (3m  -  1)  equations 


b,  -  b.  -  x.  - 
lii 


.  -  x.  +  y  (h. .  -  h, .)  p.  +  y 

i  1  lj  lj  Kj  ,L, 

j-1  J  J  J  j-1 


I,  '•  hij  *  V  “j  1  1  ‘ 2 . ”  • 


....  (21) 


e  =  r.  +  p. 

l  *i 


i  =  1 , . .  . ,m  ,  (22) 


s  .  +  q .  ; 

l  n 


i  ■  1 . m  .  (23) 


Vie  have  to  maximise  the  object  function  -y ,  given  by 


-  b  -  <-  y)  ♦  *!  ♦  l  Pi  +  I  -  hM  <Jj  (24) 


subject  to  the  constraints 

x^S-'O,  p^i^O,  q.  >  0  ,  r  ^  >  0  ,  s  ^  >  0  ;  i  -  1 , .  .  .  ,  m  .  (25) 

In  solving  this  problem  we  may  use  the  fact  that  a  feasible  solution  is 
known,  namely  0  corresponding  to  a  zero  step. 

2.6  There  are  several  ways  of  solving  linear  programming  problems,  but  in 
this  case  we  used  the  Simplex  method.  The  tableau  is  shown  as  Table  1.  The 
procedure  will  not  be  explained  hern,  but  readers  not  acquainted  with  it  will 
find  an  elementary  treatment  in  Vajda*.  It  is  an  algorithm  which  yields, 
after  a  finite  number  of  steps,  the  values  of  all  the  variables  and  the 
maximised  object  function. 

Having  obtained  p^  and  we  calculate  c.  from 

c.  -  p.  -  q i  (26) 

and  then  by  matrix  multiplication  find  the  step  vector 


z 

— o 


JL  £ 


(27) 


The  (linearised)  minimum  value  >'^u  of  y  is  obtained  by  negating  the 
maximised  object  function  of  the  linear  programming  routine. 

Fig. 2  is  a  simplified  flowchart  for  the  program. 

2  7  We  now  calculate  the  actual  value  y  of  y  at  the  new  point  of 

true  J  1 

interest.  If  it  were  not  for  ncn-linearity  this  would  be  equal  to  the  value 
ylin  8iven  by  the  linear  programming  routine.  If  c  is  small  enough  the 
two  numbers  will  be  nearly  equal. 

By  repeated  iteration  we  obtain  a  sequence  of  pairs  of  numbers 


v  (y(l)  ,  y(1)  ) 

•start  ’  ylin  ytrue; 


(y(2) 

ylin 


y(2)  ) 
ytrue;  ’ 


no 


9 


with  the  property 


y“+1)  <yU) 

1  in  ■'true 


for  all  i. 


is  small  enough 


-  (i)  ~  (i) 


so  the  sequence  {y^  is  monotone  decreasing,  approximately .  In  order  to 
guarantee  a  truly  monotone  decreasing  sequence,  we  adopt  the  following 

procedure 

"If  +  ^  >  y^*^  ,  reject  the  point  (i  +  1)  just  obtained,  replace 

•'true  'true  r 

c  by  c/2  and  repeat  the  linear  programming  routine  about  the  point  of 
interest  i". 

It  is  intuitively  obvious  that  halving  e,  if  necessary  many  times,  will  cause 

y^*  +  ^  to  approach  and  since  cannot  exceed  y^  ,  we  will 

'  true  'lin  'lin  'true’ 


eventually  arrive  at  a  value  of  which  is  not  greater  than  y'i/ 

'  'true  K  'true 

^  ( i ) 

Thus  we  obtain  a  sequence  ^ytrue^  which  is  genuinely  monotone  non-increasing, 

by  this  device  of  variable  step  length.  The  sequence  is  bounded  below  (by  0) 
and  therefore  it  converges. 

In  practice  it  was  found  thac  the  sequence  converged  quite  rapidly 
(typically  20  iterations)  until  the  values  of  v  were  constant  apart  from 

rounding  errors.  The  values  of  e  did  not  approach  zero.  In  the  tinai  steady 
state  all  the  gradients  grad  y^  are  zero,  and  all  the  y^s  equal. 

2.8  The  variable-t  device  was  also  used  to  speed  up  the  convergence,  by 
ouhling  c  when  we  appeared  to  be  f.'r  from  a  final  solution.  The  actual 
rule  adopted  was  as  follows 

,ic  -(i)  -(i+1)  -(i+2)  ,  •  i  r  ' 

Suppose  y  ,  y^  ,  yt  are  three  successive  values  of  v. 

1  'true  true  'true 

Thun  if  the  difference  between  the  last  pair  is  more  than  half  the  difference 
between  the  first  pair,  i.e. 


y(i+2>  (^>  _v(i>  ) 

'true  'true  1  'true  ‘true 


replace  r  by  2l  on  the  next  iteration”.  This  rule  allows  us  to  start 
■'iJi  a  very  snail  value  of  ,  say  10  .  The  computer  will  then  keep 


■iitliilililiialiliiiUduhLii. 
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doubling  until  it  approaches  the  final  value,  or  until  non-linearity  effects 
cause  the  's -halving'  rule  previously  described  to  come  into  play. 

No  arrangements  were  made  in  the  program  to  halt  it  when  convergence  was 
obtained.  Successive  values  of  ytrue  where  printed  and  the  program  was 
interrupted  by  the  operator  when  they  became  steady.  The  results  were  then 
obtained  via  common  storage  by  entering  an  auxiliary  program. 

3  APPLICATION  TO  2D  ARRAY 

3.1  Fig. 3  shows  the  array  consisting  of  N  omnidirectional  elements,  which 
we  label  1,  and  in  which  element  j  has  coordinates  (x  ^ ,  y^). 

The  ai ray  may  be  regarded  as  lying  in  one  plane  (:cy) .  We  shall  consider  the 
polar  diagram  for  directions  lying  only  in  this  plane. 

The  'array  function'  f (A) ,  which  is  a  complex  amplitude,  is  given  by 


where  A  »  angle  measured  from  the  x-axis 
k  ■  2 " /wavelength 

w  . . ,w  are  the  fcomglex)  weights  associated  with  the  elements. 
The  array  power  function  will  be  derined  to  be 


P(A)  -  |f (A) |2 


(32) 


3.2  The  problem  may  now  be  stated  as  follows: 

Given  N,  k,  yN,...,yN,  choose  wN,...,w;<  so  that 

the  side-lobe  level  is  minimised,  subject  to  the  condition  that  the  beamwidth 

has  a  prescribed  value  2A  . 

o 

3.3  The  restrictions  on  beamwidth  will  be  taken  to  mean  that  the  following 

equations  hold 

F(0)  =  1 

P(A  )  =  J  ,  P(“  A  )  - 

o  o 


J 


(33) 
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Thus  2A  is  the  -3  dB  beamwidth  (Fig. 4). 
o 

The  side-lobes  will  be  defined  as  all  the  maxima  of  P(A)  that  do  not 

lie  in  the  range  -A  to  A  ;  i.e. 

•  o  o 


A 


A. 

l 


where  {A^}  are  the  solutions  of 


P ' (A)  -  0  1 

P"  (A)  <  0 
Ao  <  |A|  <  n 


i  “  1 . m 


(34) 


are 


The  actual  values  of  the  side-lobes,  the  {y^}  of  the  previous  section. 


P(A-) 


(35) 


3.4  The  equation  (33)  may  be  expressed  in  terms  of  the  array  function  as 


f(0)  -  1 

f (A  )  *  2  ^  exp  (i  e  )  ► 

c  1 

f(-  A  )  ■  2  exp  (i  e^)  „ 


(36) 


by  introducing  the  (as  yet  unknown)  phases  e,  and  e (No  phase  need  be 

introduced  in  the  first  equation.)  We  then  have  three  equations  coupling  the 

N  variables  Wj,...,w^.  Provided  that  N  is  at  least  3,  we  can  solve 

these  equations  (in  general)  for  any  3  of  {w^ }  in  terms  of 

e, ,  e  ,  {x.},  {y.},  A  ,  k  and  the  remaining  {w.}.  We  choose  to  solve  for 
i  i  i  l  o  j 

w, ,  w,  and  w  .  There  are  then  N-3  'free'  variables  w.  , . . . ,w  •  or, 

1  L  J  *4  N 

rather,  since  {w^}  are  complex,  there  are  2N-6  free  real  variables  at  our 
disposal,  to  which  must  be  added  the  e^  and  giving  2N-4  free  variables. 

Thus  we  set  n  »  2N-4,  and  z^,  z^,...,z  will  correspond  to 


12 


130 


Cj,  e2,  Re(w4),  In(v4) .  Re(wN)  ,  Im(wN) 


The  equations  which  express  ,  y  in  terms  of  the  free  variables 

are  cumbersome,  but  may  be  found  in  Appendix  A. 

3.3  Having  forced  the  array  function  to  satisfy  the  mainbeam  conditions  by 
solving  the  equations  for  w^,  w0,  ,  we  now  need  to  locate  the  maxima.  The 

met hod  used  was  to  compute  P (A)  for  every  A  at  suitable  intervals 
Isay  10  di.'O  in  order  to  find  the  maxima  approximately,  and  then  refine  by 
solving 

f ' (A)  -  0  ,  (38) 

by  Newton's  method.  The  details  are  given  in  Appendix  B. 

All  solutions  with  | A !  ^  A^  are  then  deleted  and  the  remaining  angles 
re-ordered  so  that  P(A^)  is  the  largest  of  the  P(A^). 

3.6  The  procedures  described  in  section  2  call  for  the  gradients 


grad  P(AJ  , 


i  -  1,. ..  ,m 


taker,  with  respect  to  the  2N-4  dimensional  vector  z.  Now 


3P(A.) 

grad  P  (A . )  =  grad  ?(A.)  +  — — —  grad  A. 

’  A .  1 


where  grad^  denotes  the  gradient  calculated  as  if  A^  did  not  depend  on 
£.  Fortunately,  since  A^  is  a  maximum. 


■P  (A. ) 

-tst  '  » 

1 


grad  P ( A . )  =  grad  P(A.)  . 

X  Ail 

1 


and  t>n 
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The  actual  components  will  be 


[ grad  P(A.)] ^ 

[  grad  P(.A.)]  ^+1 
[grad  P(Ai)l2j42 


Dl7|P(Ai)l 

J 

3 (Re (Wj+2)J 
_ 3 _ 


(  P(A.)J  ; 

(  P(A.)]  ; 


j  -  1,2 

j  =  1.....N-3 

j  =  1 . N-3 


(42) 


(43) 


(44) 


When  we  carry  out  the  differentiations,  remembering  that  the  free  variables 
affect  P(A^)  not  only  explicitly  but  also  via  the  variables  w^,  w^,  w^ ,  the 
resulting  expressions  are  rather  lengthy,  so  these  are  relegated  to  Appendix  C, 
(in  which  the  components  of  grad  P(A^)  are  written 


di*  d2’  81’  g2 . 82N-6}  ’ 

The  numerical  computation  is  not  as  long  as  the  expressions  might  suggest  since 
the  computer  can  use  much  of  previously  stored  information. 

A  NUMERICAL  EXAMPLE 

A.l  The  example  uses  six  elements,  placed  regularly  around  a  circle  of  radius 
0.25  wavelength  (Fig. 5).  The  main  beam  is  to  be  mid-va;.  between  two  elements. 

N  =  6 

x  =  x,  *  0.216506 

1  o 

x2  =  x^  *  0 

x„  =  x.  ■  -0.216506 

3  A 

yx  =  y3  =  o-i25 

y2  =  o.25 

yA  “  y6  =  “0.125 

y5  =  -0.25 

k  * 


2  7T 
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-  The  most  natural  way  to  obtain  a  beam  in  the  x-axis  direction  is  to  apply 
phases  to  bring  the  elements  in-phase  in  this  direction;  and  then  use  equal 
weighting  an-pl  i  tudes  .  This  gjves 

wj  =  tr  exP  (-  i  k  Xj)  ;  i  =  1 , . . . ,N  (45) 

or 


w 

1 

"  W6  ’ 

0.034816  - 

0.162990  i 

w? 

m  u  is 

5 

0.166667 

w3 

a  U  a 

4 

0.034816  + 

0.162990  i 

The  polar  pattern  of  this  array  is  plotted  in  Fig. 6.  The  beamwidth  is 
84  deg.  There  are  two  side-lobes,  located  at  ±162.3  deg  with  level  -11.15  dB. 

This  pattern  should  be  compared  with  the  later  results  obtained  by  the 
side-lobe  reduction  program. 

‘•3  The  high  side-lobe  levels  under  natural  phasing  makes  this  array  a 
suitable  subject  tor  the  program,  provided  we  do  not  demand  beamwidths  much  less 
than  sh  deg.  The  actual  values  of  ?A  used  were  55  deg  to  90  deg  in  steps  of 

5  deg. 

The  initial  values  of  the  free  variables  were,  at  first,  chosen  to  be 

e  =  0,  e„  =  0 

1  L 

w  -  0.0348162  +  0.162990  i 

4 

w^  =  0.166667 

v  =  0.0348162  -  0.162990  i  , 

the  iast  three  being  taken  from  the  natural  weighting.  It  was  later  found  that 

t  h  >  '  happened  to  be  a  rather  unfavourable  starting  point,  and  machine  time 

could  ho  saver  by  starting  from  the  values  of  e.  to  w,  that  constituted  the 

i  6 

final  value e  for  anothe:  A  case. 

o 

hi,'.7  illustrates  tie  convergence  of  the  reaximum  cide-lobe  as  a  plot 
against  iteration  number.  The  step  length  is  also  shown  in  Fig. 8.  For  this 
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-4 

example,  2Ao  •  85  deg.  The  final  power  value  was  4.40172  x  10  or 

-33.56  dB  in  this  case,  showing  an  improvement  of  over  20  dB  compared  with  the 

natural  weighting  scheme  for  approximately  the  same  beamwidth. 

4.4  The  results  for  the  eight  values  of  Aq  are  given  in  Tables  2  to  9. 

For  each  element,  the  tables  give 

(i)  the  x  and  y  coordinates  in  wavelengths 

(ii)  the  weights  w^  in  real  and  imaginary  parts 

(iii)  the  weights  in  polar  form,  that  is,  amplitude  and  phase,  with  the 
phase  given  in  radians  and  in  degrees. 

The  tables  also  contain  the  results  of  an  independent  program  which 
computed  the  -3  dB  points  and  listed  the  side-lobes. 

The  actual  polar  diagrams  corresponding  to  these  results  are  plotted  in 
Figs. 9  to  16. 

It  may  be  noted  that  the  -3  dB  points  of  these  curves  occur  at  the 
required  angles.  Further,  all  the  side-lobes  are  at  the  same  level.  This 
common  level  depends  on  the  beamwidth;  the  larger  the  beamwidth  we  can  allow, 
the  lower  the  side-lobe  level .  The  trade-off  between  beamwidth  and  side-lobe 
level  is  illustrated  in  Fig. 17,  for  this  particular  array.  The  point  for  the 
original  phasing  is  also  plotted  on  this  figure,  and  it  is  about  21  dB  above 
the  curve.  As  the  beamwidth  increases,  the  side-lobe  level  drops  rapidly. 

As  the  beamwidth  decreases,  the  level  increases  as  if  to  approach  0  dB  at  about 
40°.  No  solutions  have  been  found  which  give  reasonable  patterns  for  beam- 
widths  less  than  50°,  which  suggests  that  'supergain'  weightings  do  not  exist 
for  this  particular  array. 

4.5  The  tolerances  for  these  arrays  are  also  of  interest.  Naturai.lv , 
reducing  the  side-lobe  level  makes  the  pattern  more  sensitive  to  phasing  and 
other  errors,  as  it  is  necessary  to  compare  tolerances  in  a  wav  which  is  not 
masked  by  this  effect.  In  this  Report  we  express  the  tolerance  a.  the  rms 
phase  error  which  applied  (independently)  to  all  elements,  leads  to  a  variance 
of  the  complex  array  function  equal  to  0.001.  (Rough jy  speaking,  this  means 
that  the  'noise'  in  the  pattern  is  30  dB  down.)  The  value  can  readily  be  shown 
to  be 
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If (0)1 


0.001 


N 


I  I 

J-1 


Wi' 


180  , 

— —  degrees 


(46) 


This  tolerance  is  plotted  in  Fig. 18  against  beamwidth.  On  this  curve  is  also 
shown  the  tolerance  for  natuial  phasing.  Note  that  the  tolerance  is  maximum  at 
around  84°,  the  original  beamwidth,  and  becomes  less  at  narrower  beamwidth, 
being  about  half  the  original  value  at  55°.  For  beamwidths  around  70°  to  90° 
the  tolerance  is  only  slightly  less  than  the  original  tolerance  with  natural 
phasing  (over  80%).  The  fact  that  the  side-lobe  levels  are  shown  going  down 
to  some  -40  dB  does  not  of  course  mean  that  they  could  necessarily  be  achieved 
in  practice.  It  would  be  necessary  to  consider,  for  example,  mutual  coupling, 
inter-element  screening,  departures  from  omnidirectionality,  as  well  as  the 
phase  and  amplitude  tolerances.  However  these  effects  apply  to  any  array,  and 
it  is  outside  the  scope  of  this  Report  to  include  them. 

4.6  The  weights  and  phases  for  these  optimised  arrays  have,  in  this  example, 

rather  curious  properties.  In  the  first  place,  they  are  not  symmetric  about 
the  beam  axis,  that  is 


6’ 


w2*w5. 


It  follows  that  each  solution  is  one  of  a  pair,  the  other  being  obtained  by 
reflection  in  the  x-axis.  Secondly,  the  array  has  conjugate  complex  symmetry 
about  the  y-axis,  that  is, 


‘3’ 


and  w.  and  w.  are  real. 


Probably  these  properties  are  due  to  the  geometry  of  this  example,  but  they  show 
chat  these  optimised  arrays  may  be  quite  different  from  the  arrays  one  would 
normally  consider. 

5  SUMMARY 

5.1  We  have  described  a  method  of  adjusting  the  complex  vTeights  in  a  two- 
dimensional  array  in  order  to  optimise  its  pattern.  The  optimisation  consists 
of  minimising  the  side-lobe  level  (that  is,  the  largest  side-lobe)  while  holding 
the  -3  dB  beamwidth  at  any  value  we  choose. 
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The  procedure  is  a  gradient  method,  modified  to  deal  with  the  discontinu¬ 
ities  in  the  gradient  vector.  At  each  iteration  we  choose  the  next  step  using  a 
linear  programming  routine.  The  step  length  is  controlled  by  a  parameter  which 
is  varied  automatically  so  as  to  minimise  the  number  of  steps  needed. 

T.,;  method  is  a  general  one  in  that  it  can  be  applied  to  any  problem 
requiring  the  minimising  of  the  largest  of  several  non-linear  functions. 

5.2  As  an  example,  we  considered  a  six-element  array  shaped  like  a  regular 
hexagon,  with  the  main  beam  mid-way  between  two  elements.  Starting  from  the 
natural  phasing  scheme,  the  program  reduced  the  side-lobe  level  from 
-11  to  -32  dB  for  unchanged  beamwidth.  The  solutions  indicated  weights  and 
phases  strikingly  different  from  the  natural  phasing  scheme. 

Figs. 9  to  16  show  the  polar  diagrams  for  this  array  for  various  beamwidths. 
The  minimum  side-lobe  level  decreases  quickly  as  the  beamwidth  is  allowed  to 
increase,  and  Fig. 17  shows  the  relation  for  this  particular  example.  Such  a 
curve  may  be  used  to  select  a  beamwidth  when  the.  side-lobe  level  is  prescribed. 

The  tolerance,  expressed  as  the  rms  phase  error  which  would  produce 
pattern  noise  at  -30  dB,  were  only  slightly  less  than  those  for  the  original 
phasing  (typically  80%),  so  there  is  no  objection  to  the  optimised  arrays  with 
regard  to  tolerance. 
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Appendix  A 

A.l  We  write  the  equations  in  matrix  form 


B  C] 


[:] 


(A-l) 


where  [  B  C)  is  the  6  *  2N  matrix  of  coefficients 


cos  b 
sin  b 
cos  b 

sin  b 
cos  b 
sin  b. 


sin  b 
cos  b 
sin  b 
cos  b 
sin  b 
cos  b 


cos  bj 
sin  b 
cos  b 
sin  b 
cos  b2 
sin  b. 


sin  b° 
cos  b2 


sin  b2 
.  + 


cos  b2 
sin  b2 
cob  b^ 


(A-2) 


in  which 


b° 

J 

.  + 


b. 

J 


k  x. 


k(x.  cos  A  +  y.  sin  A  ) 
J  0  J  o 

k(x.  cos  A  -  y.  sin  A  ) 
j  o  i  o 


(A-3) 

(A-4) 

(A-5) 


is  the  2N  column  vector  whose  transpose  is 


[ Re(Wj )  ,  Ira(w^) ,  Re(w„),  Im(w2) ,  ...J 


(A-6) 


E  is  the  six-vector  whose  transpose  is 


I  1  *  °>  er  (.11*2)  sin  e^  (1//2)  cos  e2>  (1//2)  sin  e2]  .  (A-7) 
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Tlie  p.i i  l  il  i on i it)’,  id  such  that  H  is  6  '  6,  C  is  6  *  (2N-6) ,  w  is 
t>  ■  1  ,  i s  (2N-b)  -  I.  This  allows  us  to  write 


A. 2  Tii.  ro  exists  the  possibility  that  b  will  be  singular,  in  which  case  the 

method  will  f.ii!.  This  will  happen  if  for  example,  A  "0.  However  for 

o 

rt-asoiKiMi-  Jit.i  no  trouble  hcs  been  encountered. 
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Appendix  B 


B.l  The  first  stage  in  locating  the  maxima  begins  by  computing  the  array 

V  ’  PM+1 

where 


p.  -  P(n(2./M  -  1)) 

PM+1  ’  P1 


i  =  1 , .  .  .  ,M 


(B-l) 


for  some  convenient  integer  M.  M  lust  be  chosen  large  enough  so  that  the 
interval  2n/M  is  less  than  the  interval  between  maxima,  and  is  found  by 
experiment.  We  then  examine  the  list  and  find  all  j  such  that 


< 


'i+1 


(B-2) 


giving  as  a  first  approximation 


A  -  r,  (2i/M  -  1)  .  (B-3) 

B.2  The  exact  value  of  A  at  the  maximum  is  then  found  by  applying  Newton's 
metnod  to  solve 


P ' (A)  -  0 


(B-4) 


If  A  is  an  approximate  solution,  the  next  approximation  A'  is 


A'  -  A  -  P’(A)/P"(A) 


(B-5) 


The  iteration  is  stopped  when 
B.3  The  formulae  for  P(A), 


A'  -  A|  <  10  4  and  A’ 
and  P"(A) 


taken  as  the  solution. 


P '  (A) 


will  now  be  derived. 
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Writing  F^  and  F^  for  the  real  and  imaginary  parts  of  f (A) ,  i 


V  . 


j-1 


h3,  Fj 


N 

k  b‘ 


where 


b.  ■  Re(v.)  cos  b  -  Im(w.)  sin  b 
3  j  o  j  o 


b,  »  ke(w.)  sin  b  +  Im(w.)  cos  b 
4  j  o  j  o 


k  x .  cos  A  ♦  k  y .  sin  A 
J  J 


then 


P(A) 


2  2 

F  +  F, 
o  1 


I  The  b.s  depend  on  j,  but  the  notation  ignores  this  for  simplicity.) 
The  first  derivative  P'(A)  is 


P'(A)  -  2< 


dF  dF,  ' 

F  — -  +  F  — - 
o  dA  1  dA 


where 


’"oVVi1 


db. 


h  dA 
J 


db 

l  -  b,  ° 

j 


4  dA 


v 


-  Vi 


(H-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 


(B-ll) 


(B-12) 


130 


Appendix  B 


23 


and 


^A 

dA 

=  l  b3  bx  (B-l 3) 

j 


where 


k  x .  sin  A  *  k  y.  cos  A 
3  j 


(B-l A) 


The  second  derivative  F"(A)  is 


and  similarly 


f  dF  dF  dF  dF 

P"(A)  =  2  <F  -rf-  +  F,  -r—  +  F,  -rr-  *  F,  — - 

]  o  dA  2  dA  1  dA  3  dA 


or 


where 


2  2 

F  F.  +  F  F  +  F  +  F- 

o  A  15  2  3 


F,  ■  ?5T<-b4bi> 


CB-15) 


db  db 

bA  dA~  ~  bl  dA~ 


=  T  [  -  b  (-  b  )  -  b  (b  b  )] 
t  4  O  1  J  1 


•  »o  -  b3  h? 


(B— 1 6) 


•  (-  b3  bo  -  b4  bj)  . 


(B-17) 
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C.l  If  i  is  any  variable 


where 


-p<V  ' 


n  r  ... 

F  =  )  Re(w.)  cos  b'^  -  Im(w.)  sin  b 

o  l-*  1  0  j 

j=i  L 


N 

\  *  l  Re<w; 
j=i  L  J 


)  sin  +  Im(w.)  cos  b 

°  J 


?■] 

?] 


(C-l) 


(C-2) 


(C-3) 


and 


f :  \ 

b',J  =  k  x.  cos  A  +  k  y.  sin  A 
o  j  'j 


3F  3F, 

r  =  0  r  _  1 
Go  ~  *  G1  ~ 


(C-A) 


(C-5) 


C.2  To  evaluate  d^,  take  a  to  be  e^.  Then 


.)  cos  -  Im(w.)  cos 

'  o  j  o 


(C-6) 


,  n  r 

G  _  =  r—  y  Re  (w. 

0  Seij=iL  J 

Now  Wj,  w 2,  w^  are  functions  of  e^  but  the  remaining  weights  are  not 


Hence 


3 

g  =  y 
0  j=i 


t— ■  (  Re(w.)l  cos  b^  -  (  Im(w.)]  sin  b^ 

Je!  J  o  3el  J  0 


.  (C-7) 


However,  Re (w^ )  is  the  (2 j  — 1 ) th  component  of  W,  and  is  equal  to 


6 

I  A 

k=l 


2j-l,kEk  (DZ)2j-l 


(C-8) 


while  Im(w.)  is 
J 


<K»«j  ■ 


(C-9) 
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Differentiating  w.r.t.  e^  gives 


T —  Re(w. ) 

j 


6  aE. 

y  a  — 

2J-1»k  3e! 


(C-10) 


— —  Im(w . ) 


3e 


1 


6  ?E, 

y  a  .  _ — 

k-T  ,k  dei 


(C-ll) 


But 


3e. 


0  unless  k  = 


3  or  4,  and  - — 
3e, 


E4>  —  -  E3  .  (C-12) 


Hence  we  get 


—  Re (w.)  -  V 

3e.  j  l,2j-l 


- —  lm(w.)  *  V,  „ . 

J  i,2j 


(C— 13) 

(C— 14) 


where 


l,k 


\,3  E4  +  \,4  E3  ’  kslto6 


(C-15) 


Similarly,  we  write 


Re (w . )  ■  V 


1 


2 , 2  j — 1 


(C— 16) 


Re  (w . )  «  V 

3e2  j  2 , 2 j  — 1 


where 


(C— 1 7 ) 


2,k 


\,5  E6  +  \,6  V  k  “  1  t0  6 


(C— 18) 


By  substitution 


j,  [V’M 


cos  -  V.  ..  sin  b^ 

o  k,2j  o 


(C— 1 9) 
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and,  similarly. 


Vk,2j-1  b< 

j"l  L 


-  V,  0 ,  cos 
k,2j  o 


(C-20) 


whore  k  -  1  for  differentiation  w.r.t.  e^,  and  2  for  ej. 

C.3  To  find  we  take  a  to  be  z^.  In  differentiating  F^  and  we 

have  to  remember  that  ,  Wj ,  depend  on  z^  via  the  matrix  D,  and  also 
that  one  of  w^,  ...  ,  W2N_^  depends  directly  on  z^. 

If  k  is  odd,  say  2 i  —  7  where  4  <  i  <  N,  then 


’  "  I  I  “  Do-  i  v  cos  ♦  D, .  .  sin  b^  I  ♦  cos  b^ 

°  j-1  L  2-1_1’k  °  2J‘k  0  J  0 

1.  *  y  I  -  D,.  .  .  sin  b  ^  ^  -  D, .  .  cos  b^  1 

1  j-i  L  2j-1*k  °  2j*k  °  j 


h^)  4.  V  (i) 

,  cos  b  +  sin  b 
k  o  o 


If  k  is  even,  equal  to  2i  -  6,  then 


1  »  [  f-  D».  .  .  cos  b^  +  D0 .  ,  sin  b^M  - 

o  2 J ~  1  »k  o  2j  ,k  oj 

1.  J  I  -  D,.  .  .  sin  b ^  -  D0 .  .  cos  b^M  + 

1  2j-l ,k  o  2j  ,k  o  J 


cos  b 


In  either  rase,  the  gradient  component  is 


8k  "  2(JoVFlV  ’ 


(C-21) 


(C— 22) 


(C-23) 


(C-24) 


(C-25) 


CM 

cO 

JZ 

\ 

4: 

i 

x: 

i 

CM 

ro 

E 

43 

43 

1 

1 

“  * 

•  •  1 

u* 

•  #  4 

•  • 

1~ 

- 

.... 

* 

43 

1 

H 

_ 

43 

X) 

43 

B  - 


max 


Table  2 


BEAMWIDTH  55  DEG 


N-  G  MAX. SI."  .  140971  (-8.50869  DB) 

El-  9.07103E-03  E2— 8.73477E-03  RAD 


1 

2 

3 

4 

X 

.216506 

0 

-.216506 

-.216506 

Y 

.125 

.25 

.125 

-.125 

WT-R 

-.28005 

.570239 

-.280266 

3.53709E-02 

wt-i 

3.08759E-02 

-2.37417E-04 

-3.03791E-02 

.190323 

WT-AMP 

.281747 

.570239 

.281908 

.193582 

PH-RAD 

3.03178 

-4 . 16347E-04 

-3.03362 

1.38705 

PH-DEC 

173.708 

-2 . 38549E-02 

-173.814 

79.4719 

5 

6 

X 

0 

.216506 

V 

-.25 

-.125 

WT-R 

.219784 

3 . 52782E-02 

WT-I 

-3.87944E-06 

-.190247 

VT-AMP 

.219784 

.19349 

PH-RAD 

-1 . 76512E-05 

-1.38744 

PH-DEG 

-1.01134E-03 

-79.4948 

1/2  PHR  -27.5 

27.5 

DEG 

0 

DB 

@ 

0 

DEG 

-8.51 

DB 

101.3 

DEG 

-8.51 

DB 

<? 

180 

DEG 

-8.51 

DB 

<?  -101.3 

DEG 

DIR  4 

.31995 

( 

6.35 

DB) 

-30DB 

TOLS: 

POS 

6.45862E- 

-03 

WL 

WT 

.352439 

DB 

PH 

2.3251 

DEG 

BEAMWIDTH  60  DEG 


N-  6  MAX.SL-  6 . 78598E-02  (-11.6839  DB) 
El-  4.96462E-03  E2— 5.02428E-03  RAD 


1 

2 

3 

4 

X 

.216506 

0 

-.216506 

-.216506 

Y 

.  125 

.25 

.125 

-.125 

WT-R 

-.202399 

.474223 

-.20235 

5 . 35152F.-02 

WT-I 

-2.20666E-02 

5. 18891E-05 

2 . 1961 7E-02 

.167044 

WT-AMP 

.203598 

.474223 

.203539 

.175407 

PH-RAD 

-3.03299 

1.09419E-04 

3.03348 

1 .26076 

PH  DEG 

-173.778 

6 , 26925E-03 

173.806 

72.2363 

5 

6 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.218211 

5 . 35400E-02 

WT-I 

-1.08309E-05 

-.167023 

WT-AMP 

.218211 

.175394 

PH-RAD 

-4.96350E-05 

-1.26059 

PH-DEG 

-2.84388E-03 

-72.2265 

1/2 

PWR 

-30 

30 

1  DEG 

0 

DB 

0 

0 

DEG 

-11. 

.68 

DB 

(3  - 

105.7 

DEG 

-11, 

.68 

DB 

<3 

105.7 

DEG 

-11, 

.68 

DB 

@ 

180 

DEG 

DIR 

4. 

86244 

( 

6.87 

DB) 

-30DB  TOLS: 

POS 

7 

.79463E- 

-03 

WL 

WT 

• 

425343 

DB 

PH 

2 

.80607 

DEG 

30 
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Table  4 


BEAMWIDTH  65  DEG 


N-  6  MAX.SL-  3.11699E-02  (-15.0626  DB) 

El-  2.90223E-03  F.2— 2 . 88797F.-03  RAD 


1 

2 

3 

X 

.216506 

0 

-.216506 

Y 

.125 

.25 

.125 

WT-R 

-.127046 

.391772 

-.127047 

WT-I 

-6 . 50121E-02 

-3.34624E-06 

6. 50226E-02 

WT-AMP 

.142714 

.391772 

.14272 

PH-RAD 

-2.66861 

-8. 54131E-06 

2.66855 

PH-DEC 

-152.9 

-4.89381E-04 

152.897 

5 

6 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.226542 

5.51492E-02 

WT-I 

-2 . 78361E-06 

-.145485 

WT-AMP 

.226542 

.155587 

PH-RAD 

-1.22874E-05 

-1.20846 

PH-DEG 

-7.04016E-04 

-6','.  2397 

1/2 

PWR 

-32.5 

32.5 

DEG 

0 

DB 

0 

DEG 

-15, 

.06 

DB 

180 

DEG 

-15, 

,06 

DB 

0  ■ 

-111.1 

DEG 

-15, 

,06 

DB 

0 

111.1 

DEG 

DIR 

4 .' 

98824 

( 

6, 

.98 

DB) 

-30DB  TOLS: 

POS 

9 

.  28274E- 

-03 

WL 

WT 

»  . 

506547 

DB 

PH 

3 

.34178 

DEG 

4 

-.216506 

-.125 

5.51484E-02 

.145494 

.155595 

1.20849 

69.2412 


BEAMWIDTH  70  DEG 


N-  6  MAX.SL-  1  33303E-02  (-18.7514  DB) 

El— 4.<*7159E-04  E2-  3.25132E-04  RAD 


1 

2 

3 

4 

X 

.2165(06 

0 

-.216506 

-.216506 

Y 

.125 

.25 

.125 

-.125 

WT-R 

-6.98724E-02 

.334035 

-6.98346E-02 

6 . 1 2264E-02 

NT- 1 

-9 . 16232E-02 

6.03756E-05 

9 . 14701E-02 

. 134584 

WT-AMP 

.115226 

.334035 

.115081 

.147856 

PH- RAD 

-2.22231 

1 . 80746E-04 

2.22286 

1.14385 

PH-DEG 

-127.329 

1.03560E-02 

127.361 

65.5378 

5 

6 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.227295 

6 . 12369E-02 

WT-I 

-4.09895E-06 

-.134572 

WT-AMP 

.227295 

.14785 

PH-RAD 

-1.80336E-05 

-1.14375 

PH-DEG 

-1.03325E-03 

-65.5321 

1/2 

PWR 

-35 

35 

'  DEG 

0 

DB 

0 

0 

DEG 

-18. 

75 

DB 

13 

117.1 

DEG 

-18. 

75 

DB 

0  - 

■117.1 

DEG 

-18. 

75 

DB 

0 

180 

DEG 

DIR 

4.81854 

( 

6.83 

DB) 

-30DB  TOLS: 

POS  1.04158E-02  WL 

WT  .568377  DB 

PH  3.74968  DEG 
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Table  6 

BEAMWIDTH  75  DEG 


N-  6  MAX. SL"  5  13755E-03  (-22.8924  DB) 

El-  4. 13922E-04  E2— 4 . 20647E-04  RAD 


1 

2 

3 

4 

X 

.216506 

0 

-.216506 

-.216506 

Y 

.125 

.25 

.125 

-.125 

vrr-R 

-1 . 75293E-02 

.288104 

-1.75279E-02 

6. 21181E-02 

wr-i 

-.110457 

3.13678E-06 

.110447 

.126866 

WT-AMP 

.111839 

.288104 

.111829 

.141257 

Plt-RAD 

-1.72818 

1.O8877E-05 

1.72818 

1.11547 

Pll-DEC 

-99.0175 

6. 23817E-04 

99.0175 

63.912 

5 

6 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.229104 

6. 21188E-02 

WT-I 

-8 . 76854E-07 

-.126864 

WT-AMP 

.229104 

.141256 

PH- RAD 

-3.82732E-06 

-1.11546 

PH- DEG 

-2.19289E-04 

-63.9114 

1/2 

PWR 

-37.5 

37.5 

DEG 

0 

DB 

0 

DEG 

-22. 

.89 

DB 

@  - 

-123.7 

DEG 

-22, 

.89 

DB 

123.7 

DEG 

-22. 

.89 

DB 

(? 

180 

DEG 

DIR 

4. 

56076 

( 

6.59 

DB) 

-30DB  TOLS : 

PCS 

1 

.12424E- 

-02 

WL 

WT 

613481 

DB 

PH 

4 

.04724 

DEG 
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Table  7 

BEAMWIDTH  80  DEG 


N“  6  MAX. SL"  1 . 69646E-03  (-27.7045  DB) 

El— 1.84984E-03  E2-  1.85059E-03  RAD 


1 

2 

3 

4 

X 

. 216506 

0 

-.216506 

-.216506 

Y 

.125 

.25 

.125 

-.125 

WT-R 

2. 73322E-02 

.253021 

2 . 7  3332E-02 

6.28844F.-02 

WT-I 

-.120433 

3 . 52960E-06 

.120423 

.124688 

WT-AKP 

.123456 

.253021 

.123486 

.139648 

PH-RAD 

-1.34763 

1 .39498E-05 

1.3476 

1.10369 

PH-DEG 

-77.2134 

7.99266E-04 

77.2119 

63.2367 

5 

6  .  .  . 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.22988 

6 . 2885bF-02 

WT-I 

-3.63328E-06 

-.124677 

VT-AMP 

.22988 

.139639 

PH-RAD 

-1.58051E-05 

-1.10364 

PH-DEG 

-9.05567E-04 

-63.2342 

1/2  PWR 

-40 

4C 

1  DEG 

0 

DB 

0 

0 

DEG 

-27.7 

DB 

0 

130.8 

DEG 

-27.7 

DB 

0  - 

130.8 

DEG 

-27.7 

DB 

0 

180 

DEG 

DIR 

4. 

28581 

( 

6.32 

DB) 

-30DB 

T0L2 : 

POS 

1 

. 16584E- 

■02 

WL 

WT 

• 

636183 

DB 

PH 

4 

.19701 

DEG 
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Table  8 

BEAMWIDTH  85  DEG 


N-  6  MAX .  SL=  4.40308E-04  (-33.5624  DB) 

Hi  —  3  .  17555K-06  F.2=  3.21454E-06  RAD 


.216506  -.216506 

.125  -.125 

6.71194E -02  6.28798E-02 

.125597  .12468 

.142407  .139639 

1.08001  1.10369 

61.88  63.2369 


5 

X  0 

Y  -.25 

WT-R  .229864 

UT-I  -3 . 93707E-06 

WT-AMP  .229864 

PH-RAD  -1.  .71278E-05 

PII-DEG  -9 . 81352E-04 


1/2 

PWR  -42.5 

42.5 

0 

DB 

*3 

0 

DEG 

-33. 

,  ^6  DB 

0 

138.5 

DEG 

DB 

0  - 

•138.5 

DEG 

-3 

DB 

9 

180 

DEG 

Dll 

36 

( 

6.05 

DB) 

-30DB  101.S  : 

POS 

1  .  17455F.- 

-02 

VI 

WT 

.  840936 

D5 

PH 

r^. 

m 

oc 

rj 

rvj 

DEG 

6 

.216506 

-.125 

6 . 28S03E-02 
-.124669 
.139629 
-1.10365 
-63.2346 


1  2 


X 

.216506 

0 

Y 

.  125 

.25 

WT-R 

6.71 199E-02 

.22631 

WT-I 

-.125609 

3 . 1 1032E-06 

WT-AMP 

.142417 

.22631 

PH-RAD 

-1  .08004 

1 . 37436E-05 

PH-DEG 

-61.8819 

7.87452E-04 

Table  9 


BF.AMWIDTH  90  DEG 


N°  6  MAX.SL-  1.07051E-04  (-39.7041  DB) 

El-  4.18666E-04  E2— 4 .17680E-04  RAD 


1 

2 

3 

X 

.216506 

0 

-.216506 

Y 

.125 

.25 

.125 

WT-R 

.10201 

.204806 

.102015 

vrr-i 

-.126915 

6.06649E-06 

.126907 

WT-AMP 

.162829 

.204806 

.162826 

PH-RAD 

-.893764 

2.96206E-05 

.89371 

PH-DEG 

-51.2089 

1 .69714E-03 

51.2059 

5 

6 

X 

0 

.216506 

Y 

-.25 

-.125 

WT-R 

.231238 

6. 34371E-02 

WT-I 

-3.07767E-06 

-.126081 

WT-AMP 

.231238 

.141141 

PH-RAD 

-1 . 33095E-05 

-1.10464 

PH-DEG 

-7.62580E-04 

-63.291 

1/2  PWR 

-45  45 

0 

DB  @  0 

DEG 

-39.7 

DB  @  -149.4 

DEG 

-39.71 

DB  @  149.4 

DEG 

-39.71 

DB  0  180 

DEG 

DIR  3. 

79087  (  5.79 

DB) 

-30DB  TOLS: 

POS  1 

. 15987E-02  WL 

WT 

632929  DB 

PH  4 

.17555  DEG 

4 

-.216506 

-.125 

6 . 34360E-02 

.126092 

.14115 

1.10468 

63.2934 
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SYMBOLS 

azimuthal  angle 

semi-beamwidth  to  -3  dB  points 

v,  at  beginning  of  iteration 

components  of  b^ 

vector  consisting  of  c  ,...,0^ 

scalars  expressing  Az  in  terms  of  the  gradients 

phases  at  -3  dB  points 

complex  array  function 
matrix  whose  rows  are  grad  y^ 

components  of  £ 

&.  &' 

components  of  h 


indices 
2n/wavelength 
number  of  side-lobes 
number  of  elements 
number  of  free  variables 

2 

array  power  function  |f(.)| 

auxiliary  variables  used  to  convert  problem  to  standard 
linear  programming  form 

complex  weight  for  element  i 

Cartesian  coordinates  of  element  i 

power  level  of  side- lobe  i 

largest  of  y^ 

initial  value  of  y 

linearised  y  after  iteration  i 

true  y  after  iteration  i 

vector  consisting  of  z^,  . .  .  .  ,wr 

free  variables,  identified  with  e^,  ej  and  the  real  and 
imaginary  parts  of  w^,...,wN 


1 
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SYMBOLS  (Cont'd) 

uz  increment  in  z  resulting  from  iteration 

c  small  constant  limiting  step  length 


grad 


gradient  operator 


matrix  transpose 
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No . 


REFERENCE 


Author 


Title, 


etc. 


1 


S  .  V  a  j  d  a 


Linear  programming  and  the  theory  of  games. 
Methuen  4  Co.  Ltd.,  London  (1960) 


Fig. 2 


£-10 


Z(I,:«Z 


Compute  sidelcbes  ot  Z(,)and  denote  by  b(1...,br 
y :  =  max  (bj  ) 


• - <^Fi rst  iteration 

Yes  9 


''Has  y 

decreased 
v  ?  y 


G  :  -Gradient 
matrix  at  z 


H  :->G  TRNfG) 


Fill  in  tableau  and 


enter  linear  programming 
routine 


£:-£/£ 


Provisional  new  point 

Z(l)-  Z  +  trn(g)#c 


Accept  new  point 

z.-r(0 


Fig. 2  Simplified  flow  chan  for  the  second  method 


v 


H+Zot  900 


-defining  points 


*1 


♦  ■49  xIO 

P(A.) 


mturol  phasing 


Beamwidth  fdeg^) 

tolerance  (-30dB  pattern  noise) versus  beamwidth 


